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We report on the extension of the experiments (P. M. Walmsley et at, Phys. 
Rev. Lett. 99, 265302 (2007)) on the decay of quasiclassical turbulence gen- 
erated by an impulsive spin- down from angular velocity to rest of superfluid 
^He in a cubic container at temperatures 0.15 K - 1.6 K. The density of quan- 
tized vortex lines L is measured by scattering negative ions. Following the 
spin-down, the maximal density of vortices is observed after time t ~ 10fi -1 . 
By observing the propagation of ions along the axis of the initial rotation, 
CN , the transient dynamics of the turbulence spreading from the perimeter of the 

container into its central region is investigated. Nearly homogeneous turbu- 
lence develops after time t ~ lOOf^ 1 and decays as L oc i~ 3 / 2 . The effective 
q{ . kinematic viscosity in T = limit is v = 0.003k, where k = 10~ 3 cm 2 s _1 is 

| the circulation quantum. 

PA CS numbers: 67.25.dk, 47.27.Gs, 47.32.-y, 47.37.+q. 



1. INTRODUCTION 

Turbulence in superfluid helium CEE1 j n ^ ne = q limit can be fully 
characterized by the positions and dynamics of thin cores of quantized vor- 
tices making a dense tangle. A useful parameter is the total length of these 
filaments per unit volume L, the mean inter-vortex distance hence being 
I = L~ l l 2 . Different spatial distributions of local density L(r) and local 
polarization of the tangle correspond to different types of turbulent flows. 
When turbulence is generated by flow on large (quasiclassical) lengthscales 
^> £, its dynamics is of great interest as the dissipative mechanisms only 
work at much smaller lengthscales. The energy is injected into large quasi- 



P. M. Walmsley et al. 



classical eddies which are a result of correlations in the polarization of vortex 
lines. While cascading down the lengthscales, this classical energy should 
be eventually converted into that of quantized flow around individual vor- 
tices. The latter have the following degrees of freedom: deformation of the 
individual filaments (Kelvin waves) and reconnections of pairs of neigbour- 
ing filaments (change in the tangle's topology). The ultimate mechanism 
of the dissipation of energy is believed to be phonon emission by Kelvin 
waves of very short wavelength. Vortex reconnections leave sharp kinks, 
thus exciting Kelvin waves. One of the modern challenges is to understand 
the nature of the non-linear cascades of energy down the classical (> £) and 
quantum (< I) lengthscales, and especially of its transfer from the classical 
to quantum cascade at the cross-over scale ~ t. 

Very recently ^ there has been progress in generating quasi-classical 
turbulence by spin-down of a non-axially-symmetric container to rest and 
detecting the decay of the vortex density L(t) through scattering of injected 
ions off vortex filaments. The latter allowed us to quantify the rate of the 
energy dissipation at different temperatures. For the given L, this rate can 
be different for different types of turbulence^, especially in the T = limit. 
It is believed that the quasi-classical turbulence achieved in the mentioned 
experiments is nearly homogeneous and isotropic. In this paper we provide 
further details on the properties of the turbulent tangles obtained in this 
way, discuss the appropriateness of the analysis, introduce a refined fitting 
procedure and summarize the current understanding of the nature of the 
energy cascade and dissipation. 

The assumptions for the analysis of the late-time decay L(t) are as 
follows^. On the classical side of wavelengths, the energy spectrum is of the 
Kolmogorov type 

E k = Ce 2 / 3 k- 5 /\ (1) 

where the Kolmogorov constant for classical fluids was found to be C ~ 1.5 
7 . Because of the confinement by the container of size d, this spectrum 
is truncated at short wavenumbers by k\ ~ 2n/d making d the energy- 
containing length. 

At short length scales < t the cascade becomes quantum as the dis- 
creteness of the vorticity in superfluids adds new behavior. At sufficiently 
low temperatures, it is expected to take the form of a non-linear cascade 
of Kelvin waves The most complicated regime is the transitional region 
between these two limits, the Kolmogorov and Kelvin-wave cascades. The 
question being currently debated is whether the energy stagnates at k < I" 1 
due to the poor matching in the kinetic times of the two cascades ^ or gets 
efficiently converted from classical eddies to waves along quantized vortex 
lines with the help of various reconnection processes 
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In the steady state, the energy flux e is equal to the dissipation of 
the kinetic energy. The redistribution between different lengthscales and 
eventual dissipation of the flow energy is through the motion of vortex lines. 
The flux down the classical cascade of this energy per unit mass can be 
assumed to be ^ 

e = u(kL)\ (2) 

where k 2 L 2 is an effective total mean square vorticity and the "effective kine- 
matic viscosity" v can vary depending upon temperature and type of flow. 
With regard to the physical meaning of Eq. (f2]) there are two possibilities. It 
can be simply an analog of the classical formula for the viscous dissipation 
rate in shear flow at the smallest dissipative lengths. Alternatively, it can 
result from a cascade-specific process (like "bottleneck") limiting the energy 
flux e down the lengthscales. 

Equations Q] and [2] lead to the late-time free decay ® 

L = {ZCf^K-^u- 1 ' 2 ^/ 2 . (3) 

In the rest of this paper we show that this regime is observed for quasi- 
classical turbulence in superfluid 4 He generated by various means, including 
the T = limit it is also observed for vortex tangles in superfluid 

3 He-B generated by a vibrating grid 



2. Experimental 

To detect the presence of vortices we used the fact that injected ions 
interact with a vortex. This process is characterized by the trapping (or 
"scattering") diameter a which depends on temperature, electric field and 
type of ion. Two types of charge carriers have been used: at T > 0.8 K, free 
ions (an electron self-localized in a bubble), and at T < 0.8 K small charged 
vortex rings (CVRs) with one negative ion trapped on their cores. 

For free ions, the trapping diameter a is inversely proportional to the 
electric field, a oc E . In a field E = 10 V/cm it is decreasing with temper- 
ature from a = 220 nm at T = 1.6 K to a = 4nm at T = 0.8 K. In contrast, 
CVRs have a much larger trapping diameter of order their radius (typically, 
a = 0.4-1.7 (im in our experiments) which makes them much more sensitive 
to tangles of small density L. 

After injecting a short pulse of ions, a sharp pulse of current arrives at 
the collector as shown in Fig.[TJ The interaction with vortices, that happend 
to be in the way of the propagating ions leads to the depletion of these pulses. 
The decrease in the peak height I of the current pulses due to the ballistic 
CVRs (shown in Fig. [I]), relative to the long time amplitude I(oo) when all 



P. M. Walmsley et al. 



vortices have decayed, is the measure of blocking the paths of some ions 
with the probability per unit length La. We hence determine the average 
vortex density between the injector and collector, after the turbulence has 
been decaying a time t, as 

= (ad)- 1 In (^M) . (4) 

To calibrate the value of a for different types of ions and at different tem- 
perature and electric field, we measured the attenuation of the pulses of the 
collector current in the horizontal direction when the cryostat is at continu- 
ous rotation Q2HS] arounc l the vertical axis at angular velocity (thus having 
an equilibrium density of rectilinear vortex lines L = 20/ k) as in Fig.0 




Time after start of injection (s) 



Fig. 1. Examples of current records produced by short pulses for three 
different temperatures (the temperatures, pulse durations and mean driving 
fields are indicated in the legend). The solid lines show the records without 
a vortex tangle in the ions' path, while the dashed ones represent the records 
suppressed by the vortex tangle which has been decaying a specified time 
(indicated near curves) after stopping generation. The charge carriers are 
either free ions (T = 1.60 K and T = 0.90 K) or charged vortex rings (T = 
0.15 K). The electronics time constant is 0.15 s, hence the true time of arrival 
of the ions producing the fastest peak (T = 0.90 K) cannot be resolved. 



At low temperatures, Kelvin waves of a broad range of wavelengths 
are excited, thus making the tangle length a fractal property. However, 
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Fig. 2. Dependence of pulse amplitude on the angular velocity of rotation, Q, 
in examples of calibration measurements. The temperatures, driving fields 
and trapping diameters a are indicated. The charge carriers are either free 
ions (T = 1.60 K) or charged vortex rings (T = 0.78 K and T = 0.50 K). 



the main contribution to the total length L converges quickly at scales just 
below £. Hence, a probe ion with the trapping diameter a <C £, moving 
at a speed (v > 5cm/s), much greater than the characteristic velocities 
of vortex segments (~ k/£ < 10 _1 cm/s), should sample the full length L. 
Experimental measurements with different a = 0.4-1.7 /xm (using charged 
vortex rings in a range of driving electric fields) indeed produce consistent 
values of L. 

In order not to introduce extra turbulence by the probe ions, the mea- 
surement was performed by probing each realization of turbulence only once 
- after a particular waiting time t during its free decay. Then the contam- 
inated tangle was discarded and a fresh tangle generated. The values of L, 
measured with the use of either of the two types of ions in the temperature 
window 0.7-0.8 K where they co-exist, were consistent too. 

The experimental cell had cubic geometry with sides of length d = 
4.5 cm. A schematic drawing of the cell is shown in Fig.O The relatively 
large size of the cell was important to enhance the efficiency of ion trapping 
and the time resolution of vortex dynamics, and was really instrumental in 
ensuring that the continuum limit f C d holds even when the vortex line 
density drops to as low as just L ~ 10cm~ 2 . This also helped ensure that 
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Fig. 3. A side cross-section of the Manchester cell. The distance between 
opposite electrode plates is 4.5 cm. An example of the driving electric field is 
shown by dashed lines. Such a configuration was used in all the Manchester 
experiments described in this section. It was calculated for the following 
potentials relative to the right electrode: left electrode at -90 V, side, top 
and bottom electrodes at -45 V, and the right electrode and collector at 
producing a 20 V/cm average driving field in the horizontal direction. To 
inject ions, the left tip was usually kept at between -500 to -350 V relative 
to the left grid. When ions traveling vertically across the cell (injected from 
bottom tip and detected at the top collector) were required, the potentials 
on the electrodes were rearranged as appropriate. 



the presence of the walls, which might accelerate the decay of turbulence 
within some distance ~ £, does not affect the dynamics of the turbulent tan- 
gle in the bulk of the cell. In order to probe the vortex densities along the 
axial and transverse directions, there were two independent pairs of injec- 
tors and collectors of electrons. All injectors and collectors were protected 
by electrostatic grids, enabling injecting and detecting pulses of electrons. 
The injectors were field emission tips made of 0.1mm diameter tungsten 
w j re OH The threshold for ion emission was initially ~ -100 V and -210 V 
for the bottom and left injectors respectively; however, after over two years 
of almost daily operation they changed to some -270 V and -520 V. The fact 
that the two injectors had very different threshold voltages helped investi- 
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gate the dependence (or rather lack of it) of the radius of initial charged 
vortex rings on the injector voltage. The six side plates (electrodes) that 
make up the container cube can be labeled "top" , "bottom" , "left" , "right" , 
and two "side" electrodes. The top, bottom, left and right electrodes had 
circular grids in their centers. All grids were made of square tungsten mesh 
with period 0.5 mm and wire diameter 0.020 mm, giving a geometrical trans- 
parency of 92%. The grids in the bottom, left and right plates had diameter 
10 mm and the grid in the top plate had diameter of 13 mm. The injector 
tips were positioned about 1.5-2 mm behind their grids. The two collector 
plates were placed 2.5 mm behind their grids and were typically biased at 
+10 — |-25 V relative to the grids. Further details can be found in^. 

The novel technique of generating quasiclassical turbulence, suitable for 
any temperatures down to at least 80 mK, relied on rapidly bringing a rotat- 
ing cubic-shaped container of superfluid 4 He to rest ™. This was achieved by 
mounting the experiment on the Manchester rotating cryostat USED whose 
angular velocity and acceleration can be controlled over a wide range. In our 
experiments, the range of angular velocities of initial rotation f2 was 0.05- 
1.5rad/s. In classical liquids at sufficiently high Reynolds numbers, even in 
cylindrical containers, spin-down to rest is always unstable: within a few 
radians of initial rotation upon an impulsive spin-down to rest, a nearly ho- 
mogeneous turbulence develops with the energy-containing eddies of the size 
of the container. In a cubic container, the turbulence becomes homogeneous 
much faster. One expects similar processes to occur in a superfluid liquid 
providing the process of initial multiplication of vortices does not affect the 
dynamics. As spin-down experiments always begin from already existing 
dense rectilinear vortex arrays of equilibrium density L = 2Q/k, and rapid 
randomization and multiplication of these vortices is expected due to the 
lack of axial symmetry of the container, as well as surface pinning (friction) 
near container walls, these seem to be sufficient for the superfluid to mimic 
the large classical turbulent eddies. 

Before making each measurement, the cryostat was kept at steady ro- 
tation at the required f2 for at least 300 s before decelerating to full stop, 
then waiting a time interval t and taking the data point. The deceleration 
was linear in time taking 2.5 s for f2 = 1.5rad/s and 0.1s for f2 = 0.05rad/s. 
The origin t = was chosen at the start of deceleration. 

We have also found that currents of injected ions leave vortex tangles 
behind that can also be used as an alternative means of generating turbu- 
lence. Depending on the conditions of injection, such as the temperature 
and duration of the current pulse, the tangles can have features of either 
quasi-classical (large-scale dominated) or quantum (quantum energy domi- 
nated) turbulences. We do not discuss these experiments in this paper but 
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quote the resulting values of v(T) in later sections. For more details see^. 



3. Results 
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Fig. 4. Records of the current to the top collector, injected from the bottom 
injector as a 0.1 s-long pulse at the time t after an impulsive spin-down from 
n = 1.5rad/s to rest. T = 0.15 K, E = 20V/cm. 



To illustrate what is happening near the vertical axis of the container 
at different stages of the transient after a spin-down, in Fig. 0] we show five 
records of the current from the bottom injector to the top collector arriving 
after a 0.1 s-long pulse of probe ions was fired from the bottom injector. 
Each is characteristic of a particular configuration of vortex lines near the 
rotational vertical axis of the container during the transformation from an 
array of parallel lines to a homogeneous decaying vortex tangle. One can see 
three different characteristic times (vertical dashed lines in Fig. [4]) of arrival 
of ions via different means. The first peak at ss 0.4 s (determined by the time 
constant of the current preamplifier, 0.15 s) corresponds to the ions trapped 
on the rectilinear vortex lines which can slide along those lines very quickly, 
provided these lines are continuous from the injector to collector as, for 
example, during steady rotation. Our estimates show that at SI = 1.5rad/s 
of order half of rectilinear vortices from the equilibrium array terminate on 
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the collector grid and hence do not allow trapped ions to reach the collector. 
The second peak at ~ 1 s corresponds to the coherent arrival of the ballistic 
charged vortex rings from the bottom to the top. The third broad peak at 
times ~ 3 s (but with a long tail detectable until ~ 40 s) corresponds to the 
charge trapped on the vortex tangle and drifting with the tangle. Hence, we 
have the following regimes (and curves in Fig.H]) labeled by the time t after 
the spin-down: 

1. £lt < 0, steady rotation. The nearly equal first and second peaks 
indicate that there exist an ordered array as the ions can arrive at the 
collector by either fast sliding along the rectilinear lines or propagating 
as slower CVRs between the lines. 

2. £lt = 1.5. The first peak gets enhanced three-fold while the second 
one is still there (and no third peak) - meaning more trapped ions can 
now reach the collector along the rectilinear vortex lines (apparently 
because many vortex lines no longer terminate on the grid after being 
"shaken off" by the emerging turbulent layer near the horizontal walls) 
while there is not much turbulence in the central region yet. 

3. £lt = 6. The first peak has disappeared in favor of the second one 
which got broadened - at this stage the rectilinear vortices should 
have become scrambled in the Ekman layers near the top and bottom 
walls while the ballistic charged vortex rings are still the dominant 
transport of charge. 

4. £lt = 12. Now both fast peaks have disappeared completely while the 
third broad peak has emerged - this means that a dense turbulent 
tangle has finally reached the axial region. 

5. £lt — > oo. The sharp second peak has recovered but all others van- 
ished - after the turbulence has decayed only ballistic charged vortex 
rings carry the charge, neither the rectilinear array or turbulent tangle 
contribute to the transport any more. 

In Fig.[5l the measured densities of vortex lines along the horizontal 
(transverse) direction, Lt, are shown for four different initial angular veloci- 
ties fi. During the transient, which lasts some ~ lOOfi , L t (t) goes through 
the maximum after which it decays eventually reaching the universal late- 
time form of L oc t~ 3 / 2 . For 0, > 0.5rad/s, the values of L at maximum 
were too high to be detected. The initial vortex densities at steady rotation, 
L = 2Q//C, are shown by horizontal lines. Similarly, the densities of vortex 
lines measured along the vertical (axial) direction, L a , are shown in Fig.0 
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Fig. 5. L t (t) at T = 0.15 K for four values of 0. The average driving fields 
used for f2 = 1.5rad/s: 5V/cm (o), lOV/cm (y), 20V/cm (o), 25V/cm 
(A). At other values of O the electric fields used were either 20V/cm 
(0.05rad/s) or 10 V/cm (0.5 and 0.15rad/s). The line shows the dependence 
£-3/2 _ rp^g jjoj^ontal bars indicate the initial vortex densities at steady rota- 
tion, L = 2Q/k, at ft = 1.5rad/s, 0.5rad/s, 0.15rad/s and 0.05rad/s (from 
top to bottom). 



To stress the scaling of the characteristic times with the period of initial ro- 
tation f^ 1 and the universal late-time decay oc i _3//2 , we re-plot these data 
for different f2 in Fig. [7] rescaled accordingly. 

We can see that at all the transients are basically universal. Immedi- 
ately after deceleration, L t increases, indicating the appearance of the turbu- 
lent boundary layer at the perimeter, while L a is stable at L\ ~ 20/k. Only 
at ~ 3, the latter starts to grow, signaling the destruction of the rotating 
core with vertical rectilinear vortices. After passing through a maximum at 
£lt = 8 and Vtt = 15 respectively, Lt and L a merge at fit ~ 30 and then 
become indistinguishable. This implies that from now on the tangle density 
L is distributed nearly homogeneously. Eventually, after Qt ~ 100, the de- 
cay takes its late-time form L oc i~ 3 / 2 expected for quasiclassical isotropic 
turbulence, whose energy is mainly concentrated in the largest eddies bound 
by the container size d, but homogeneous on smaller length scales. We hence 
assume that the turbulence in 4 He at this stage is nearly homogeneous and 
isotropic, and can apply Eq. (|3|) to extract the effective kinematic viscosity 
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Fig. 6. L & (t) at T = 0.15 K for four values of f2. The average driving 
field used was 20V/cm in all cases except at 0.5rad/s where both lOV/cm 
(y) and 20V/cm (A) were used. The line shows the dependence i~ 3 / 2 . 
The horizontal bars indicate the initial vortex densities at steady rotation, 
L = 2Q/k, at O = 1.5rad/s, 0.5rad/s, 0.15rad/s and 0.05rad/s (from top 
to bottom). 



v (see Discussion). 

At all temperatures the transients Lt(t) after a spin-down are, in first ap- 
proximation, universal, i. e. the timing of the maximum is the same Qt ~ 7 
and the amplitudes of the maximum are comparable. This supports our 
approach to turbulent superfluid helium at large length scales as to an invis- 
cid classical liquid, agitated at large scale and carrying an inertial cascade 
down the lengthscales, independent of temperature. On the other hand, as 
the temperature rises and the damping through the normal component in- 
creases, changes in the shape of transients might be expected. Indeed, one 
can see that the slope of L(t) after the maximum but before reaching the 
ultimate late-time decay L oc t _3//2 (i.e. for 10 < Qt < 100) is changing 
gradually with increasing temperature from being less steep than i -3 / 2 at 
T = 0.15 K to more steep than t~ 3 / 2 at T = 1.6 K (Fig. [8]). At temperatures 
around T = 0.85 K (Fig. El) it nearly matches i~ 3 / 2 , thus making an erro- 
neous determination of v possible by taking this part of the transient for the 
late-time decay L oc i -3 / 2 . Indeed, in the first publication 4 , parts of some 
transient at T = 0.8 — 1.0K for as early as £lt > 15 were occasionally used 
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Fig. 7. fr 3 / 2 L t (» (filled symbols) and fT 3/2 L a (t) (open symbols) vs. Qt 
for four values of f2 at T = 0.15 K. The straight line oc i~ 3 / 2 guides the eye 
through the late-time decay. 



to be fitted by L oc i~ 3//2 that often resulted in overestimation of the value 
for v. To rectify this, we have re-fitted the datasets used in the original pub- 
lication as well as subsequent measurements with L oc i~ 3 / 2 for spin-downs 
using the following rules: for T < 0.5 K only points for £lt > 300 were used, 
between 0.5 K and 1.0 K only points for £lt > 150 used and at T > 1.0 K 
only > 75 were used. This resulted in slight systematic reduction of the 
extracted values of v(T) at temperatures 0.8-1.2 K from those published in 
4 ; what looked as a rather steep drop in v(T) now occurs at 0.85 - 0.90 K 
and is somewhat reduced in magnitude. 

4. Discussion 

Let us begin from commenting on the significance of observation of 
L oc t -3 / 2 decay in a broad range of decay time t. For quasi-classical turbu- 
lence, at lengthscales approaching the cross-over length ~ £, below which the 
classical cascade no longer holds, there might exist an additional "equiparti- 
tioned" component of vorticity that does not contribute to the steady-state 
energy cascade Also, in any type of quantum turbulence, unlike the 
classical case, there exists an additional quantum contribution to the to- 
tal energy resulting from the quantized flow around individual filaments, 
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f(s) 

Fig. 8. Free decay of a tangle produced by an impulsive spin-down to rest ^ 
from 1.5 rad/s and 0.5 rad/s, as well as by a jet of free ions from the bottom 
injector (•) 5 , at T = 1.60 K. All tangles were probed by pulses of free ions 
in the horizontal direction. The line L oc i~ 3 / 2 corresponds to Eq. ([3]) with 
v = 0.2k. 



Eq = BL where B ~ const ^. During the free decay, both these reservoirs 
release additional energy to be dissipated. For the math leading to Eq. fl5J) 
to hold, these energy fluxes should not add to the flux e of classical energy 
in Eq. ((2J) (or should at least have the same time-dependence oc t~ 3 as the 
classical one, which is unlikely). At early stages of the decay \E q \ is always 
small compared to e, but estimates show that it might become comparable 
at very late times, especially in the T = limit. The fact that we never ob- 
served deviations from the L oc t -3 / 2 decay at late times suggests that these 
fluxes are simply not coupled to Eq. ((2J). For instance, this would be the case 
if this equation describes some L-dependent limitation for the energy flux 
approaching the classical-quantum crossover scales ~ t as the quantum flux 
\Eq\ can only be injected at shorter lengthscales < I. 

In Fig. [10] we plot the values of v(T) as extracted from the late-time 
decays L oc i -3 / 2 using Eq. We assumed that, at all temperatures, 
C = 1.5 and k\ = 2n/d, and the only temperature-dependent parameter in 
Eq. ([3j) is v{T). In principle, the interaction between the vortex tangle and 
container walls changes with temperature (at high temperatures, via mu- 
tual friction with the normal component which is viscously coupled to walls, 
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f(s) 

Fig. 9. Free decay of a tangle produced by an impulsive spin-down to rest HI 
from 1.5 rad/s, as well as by a jet of free ions from the bottom injector (•) ^, 
at T = 0.85 K. All tangles were probed by pulses of free ions in the horizontal 
direction. The ion jet data are the average of nine measurements at each 
particular time but the spin-down data show individual measurements. 



but at low temperatures, only through the non-axially-symmetric shape of 
the container, including both its cubic shape and surface roughness). Hence, 
the boundary conditions of the quasi-classical turbulent flow can change with 
temperature, thus affecting the numerical prefactor in the relation for the 
cut-off wavenumber for the largest eddies, ki = 2ir/d. Also, at low temper- 
atures, the hard core of the giant vortex might survive for longer after the 
spin-down. This would result in an anisotropic rotating turbulence whose 
energy-containing eddy (the giant vortex) might release its energy slower 
than the classical isotropic 3d turbulence: as a result the effective value of 
the constant C might decrease with decreasing temperature. The fact that 
we have shown experimentally that the vortex density L is distributed homo- 
geneously throughout the cell at late times does not guarantee that no such 
large-scale anisotropy exists. Also, there is no proof that the coarse-grained 
hydrodynamic equations have the same strength of the non-linearity as the 
classical Navier-Stokes equation; hence, the effective Kolmogorov constant 
can in general differ from the classical value (and can depend on temperature 
too). While we believe that these changes with temperature of the effective 
parameters C and k\ might be unlikely, we cannot rule them out at this 
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Fig. 10. The effective kinematic viscosity v(T) for quasiclassical turbulence 
in 4 He: v(T) infered using Eq. ([3]) from the free decay L oc i~ 3//2 of tangles 
produced by impulsive spin-down (open triangles'^), ion jet (open circles^) 
and towed grid (open squares The dashed curve is a fit to our data 
(open triangles) by theory For comparison, the values for v(T) for the 
quantum regime of turbulence in 4 He, generated by a beam of colliding 
CVRs, are shown by crosses^. 



stage. Obviously, further spin-down experiments in containers of different 
shapes as well as with the turbulence generated by other means (e.g. ro- 
tating agitators towed grid vibrating grid propagating turbulent 
front counternow ^ ion jet^) and numerical simulations of the tangle in 
a cubic cell upon a spin-down can clarify these issues. For example, the fact 
that quasi-classical tangles, generated by ion jets in a non-rotating cryostat 
^, yield the same values of L(t) oc t~ 3 / 2 (see Figs.lBHTDJ), and hence of v(T), 
down to T = 0.7 K as those after spin-down strongly suggests that slowly 
decaying background rotation surviving long time after a spin-down, if any, 
might be irrelevant. 

For comparison, in Fig. llOl the high-temperature data for quasi-classical 
turbulences generated by towed grids and counternow in channels with square 
cross-section are shown too. Generally they are in good agreement with our 
results. We also show the low-temperature values of v ~ 0.1k for quan- 
tum tangles generated by short injections of CVRs, that apparently have 
little large-scale flow. The main observation is that the infrerred values of 
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v{T) for quasi-classical tangles show a substantial decrease with cooling be- 
low T ~ 1 K, eventually reaching the temperature-independent plateau of 
v ~ 0.003k that is much smaller than that for the quantum tangles. 

At the moment there exist two different theoretical models explaining 
the dropping of quasiclassical v(T) approaching the T = limit. The "bot- 
tleneck" scenario by L'vov, Nazarenko and Rudenko ^ explains the pile- up 
of excess vorticity on the classical side of the cross-over lengthscale I as a 
result of poor matching of the kinetic times of the classical and Kelvin-wave 
cascades; this pile-up can only happen at sufficiently low temperatures when 
the cascading energy can reach the Kelvin-wave cascade. On the other hand, 
Kozik and Svistunov'm explain the enhancement of the total length L at low 
temperatures as a result of extra contribution at lengthscales shorter than 
I caused by the fractalization of the tangle after various reconnections pro- 
cesses become possible at low temperatures where Kelvin waves are hugely 
underdamped. It seems, at this stage we cannot discriminate conclusively 
between the two alternative explanations; it is possible that both mecha- 
nisms contribute. Future experiments as well as numerical and analytical 
calculations should help clarify the issue. For example, it is important to 
know the rate of reconnections as function of the net polarization of the 
tangle. We are currently conducting experiments on the dynamics of tangles 
at continuous rotation that might advance our understanding. 

To conclude, we would like to note that these are the first measurements 
on superfluid turbulence below 1 K after fifty years of research which now 
point out (i) how to create homogeneous isotropic turbulence (at least in a 
transient state) and (ii) how to measure it. The road is now open to study 
the remaining open questions in the T = limit, starting from a convincing 
experimental proof of equation [2] on which the analysis in the paper is based. 
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